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Abstract 

The reliability of the models describing the Final-State Interac- 
tions (FSI) in (e,e'p) scattering at high proton energies is an impor- 
tant issue in view of the experiments planned at CEBAF. One of the 
most popular approaches adopted, the Glauber method, involves the 
linearization of the wave equation for the ejected proton travelling 
through the residual nucleus. We have studied the consequences of 
such an assumption for the case of the ^^C(e,e'p)^^B* reaction at high 
proton momenta by comparing the results with the predictions ob- 
tained when the second-order differential equation for the proton wave 
is solved exactly for each partial wave. We find that the two methods 
give well-correlated angular distributions for momenta in the range 
1-4 GeV/c, i.e. for kinematics relevant to the transition to the eikonal 
regime. 



In the experiments planned at CEBAF concerning (e,e'p) scattering [Q], 
where the proton momentum can be larger than 1 GeV/c, a central issue is 
the reliability of models describing the Final-State Interactions (FSI) between 
the outgoing proton and the residual nucleus. For example, at moderate 
missing momenta pm (Pm = P^Q^ where p is the outgoing proton momentum 
and q is the momentum transferred to the target by the electron) an accuracy 
within 10% is required to unambigously identify exotic effects like Colour 
Transparency 0, if any. 
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The most widely used approach to the problem of FSI at high energy is 
the Glauber approximation Q], since there is a long well-established tradition 
in the analysis of proton-nucleus elastic scattering [|, 0. In particular, in 
that context it has been shown that the validity of such an approach can 
arise from a non trivial cancellation among the leading corrections to the 
lowest-order theory 0. However, the generalization to the (e,e'p) scattering 
is not straightforward because the initial proton state and, in general, the 
kinematics are quite different. 

In the Impulse Approximation, the basic ingredient of completely exclu- 
sive (e,e'p) reactions is the scattering amplitude []^, H 

J^(g) = y'drdaei«V)*(r,a)>(g,r,a)vI/„(r,a) , (1) 

where x^~'^ ^ind \E'q, describe the scattering- and bound-state wave functions 
of the nucleon knocked out from a hole with quantum numbers a, respec- 
tively. Usually, the current operator J'^ is approximated by a nonrelativistic 
expansion in powers of the inverse nucleon mass, thus introducing uncertain- 
ties which become more important with increasing energy |T^. However, 
our interest is not focussed on the comparison with experimental data, but 
on the analysis of the scattering wave x^~^- Therefore, we have considered 
the simplified picture where we retain just the longitudinal component of J'^ 
in the leading order o(l) of the nonrelativistic approximation and we neglect 
the nucleon form factor. The cross section becomes, therefore, proportional 
to 

/drdaei«V-^*(r,a)vI/„(r,a)f = , (2) 

which is traditionally identified as the "distorted" spectral density [ITT 



at the missing energy corresponding to the knockout hole a. 

In the framework of the Distorted- Wave Impulse Approximation (DWIA) 
[0, |] the scattering wave function x^~^ is solution of the Schrodinger equation 

(-|^V^ + X = i^cmX , (3) 

where m is the reduced mass of the proton in interaction with the residual 
nucleus, -Ecm is its kinetic energy in the cm system and V contains a lo- 
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cal equivalent energy- dependent optical potential effectively describing the 
residual interaction. 

Eq. can be solved for each partial wave of x^~^ up to a maximum angu- 
lar momentum I/max (p), where a convergency criterion for the partial- wave 
expansion is satisfied. The boundary condition is such that each incoming 
partial wave coincides asymptotically with the corresponding component of 
the plane wave associated to the proton momentum p. Typically, this method 
(from now on method I) has been applied to (e,e'p) scattering for proton mo- 
menta below 0.5 GeV/c and Lmax < 50 Pj. In the kinematics explored in 
this work a maximum Lmax = 120 has been used. 

At higher energies the Glauber method suggests an alternative way 
(from now on method II) of solving eq. (^ by linearizing it along the prop- 
agation axis z: 



r = z^ + b (4) 

(5) 





(6) 

where 6 describes the degrees of freedom transverse to the motion of the 
struck particle with momentum p. With this approximation eq. (|^) becomes 

The boundary condition is of incoming unitary flux of plane waves. 

Both methods I and II solve the Schrodinger equation for the nucleon 
scattering wave. Relativistic effects are correctly taken into account only in 
a proper calculation of the kinematics. In the case of the application of the 
Glauber approach to unpolarized proton-nucleus scattering, this approxima- 
tion does not seem to produce relevant consequences 0, ||, |^. Nevertheless, 
because of the previously mentioned differences, the generalization to the 



(e,e'p) should be tested P, O, p 
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The eikonal approximation is supposed to reproduce the exact solution 
of eq. (Q) for pm <C g, and in general its reliability increases with the ejectile 
energy ideally in the limit where x^~^ is expanded on an infinite number 
of partial waves. On the other hand, method I can be considered reliable only 
for nucleon energies such that the condition Lmax ^ -^target P fulfilled. 

In a previous work we have analysed the spectral density of eq. 



for the ^^C(e,e'p) Bsi/2 reaction in parallel kinematics (p along q). In 
the range 1 < p < 2 GeV/c and for Lmax = 120 we found a good correlation 
between the predictions of methods I and II, the (small) discrepancies being 
ascribed to the impossibility of taking into account in eq. the interference 
between the incoming and the reflected flux. The consequent overestimation 
of S'P with respect to method I was found to be related to kinematics and 
proportional to the absorption strength of the optical potential. 

In the present work we have extended the calculations up to p, g = 4 
GeV/c by improving the numerical precision of the FORTRAN code. We have 
considered the ^^C(e,e'p)^^B* reaction in the socalled perpendicular kinemat- 
ics, where p and q are kept constant and the angle between their directions, 7, 
is allowed to vary. The bound state in eq. (§) is a solution of the Woods- 



Saxon potential of Comfort and Karp |15| with the quantum numbers a of 
the s wave. For sake of simplicity, the contribution to V coming from the 
Coulomb potential has been neglected to avoid numerical problems related 
to the high angular momenta required. Therefore, in proper terms the results 
presented here refer to the (e,e'n) reaction. V is an optical potential of the 
simple form 



V{r) = (U + iW) 



r-R 

1 + e~ 



= {U + iW)p{r) , (8) 

with R = 1.2 X A^^'^ fm and a = 0.5 fm. The nuclear density p(r) defined in 
eq. (|]) is normalized such that p(0) = 1. 

At the nucleon energies here considered, the parameters U, W can only be 
guessed. According to the Glauber model the imaginary part should scale as 
W ^ p/10 MeV, while U/W should equal the ratio between the real and the 
imaginary parts of the average proton-nucleon forward-scattering amplitude. 



which is expected to be small in the considered kinematics |]16 . 
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In fig.l S^^2 is calculated by method I as a function of V for pm = and 
p = q = 1.4 GeV/c, where the elementary cross section for the rescattering 
of the ejectile is predominantly inelastic. The three curves correspond to 
W = 0, 50 and 100 MeV, for U ranging continuously from —50 to +50 
MeV. The {U = 0, W = 0) point corresponds to the Plane- Wave Impulse 
Approximation (PWIA) result, where any rescattering between the ejectile 
and the residual nucleus is neglected. The middle curve (W = 50 MeV) shows 
an average 40% damping with respect to the PWIA result in agreement with 
the observation of the NE18 experiment in the context of a semi-inclusive 



(e,e'p) reaction at small 1131 • It is evident that the sensitivity to both 
the sign and the magnitude of the real part of the potential is very small but 
for huge values oi U ^ W, which are forbidden by the mainly absorbitive 
character of the proton-nucleon amplitude at these kinematics. Therefore, in 
this work we will use U = 0. 

Eq. ([^ produces by definition results that only depend on V through the 
ratio V/p. If \^ is a linear function of p (as suggested by the Glauber model), a 
constant absorption is produced at any value of momentum transfer, provided 
that p = q and pm — (i.e. small 7). In fig. 2 we show that the same 
property holds, with a good approximation, also for method I. The distorted 
spectral density is given as a function of p = q for pm = and the choice 
U = Q,W = p 50/1400 MeV, which produces at p = g = 1.4 GeV/c the 
40% absorption observed in the NE18 experiment. The flat curve shows 
that this damping remains constant down to very low values of p = g (p = 
q ^ 0.4 GeV/c). However, the choice of an absorption-dominated potential 
is reasonable only in the kinematical region where the elementary proton- 
nucleon scattering amplitude is dominated by inelastic processes, i.e. for 
p > 1 GeV/c |]T6|. Only above this threshold, which is anyway relevant to 



the kinematics explored at CEBAF, our considerations can be applied (for 
applications of the Glauber approach to (e,e'p) reactions at lower energies see 
refs. [1T8| , p!9[|). Since the NE18 observation of a roughly g-independent 40% 
damping at small pm is reproduced by both methods I and II with the optical 
potential of eq. (g) and with U = 0, W = p 50/1400 MeV, we adopt this 
choice in the following also for the calculation of the angular distributions. 

In the pure Glauber theory W is not a free parameter and would result 
in a larger value. As previously mentioned, the increase of W leads to a 
larger discrepancy between the results of methods I and II, but the overall 
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agreement is not too much spoiled |T^. However, the 14^- value suggested by 
the Glauber model produces a sensibly larger absorption than observed in 
the NE18 experiment (for more details on this topic see refs. |T^, ^ p2|). 



In figs. 3 and 4 the distorted spectral density 5*1^/2 is shown as a function 
of 7 for g = p = 1.4 and 4 GeV/c, respectively. Because of the large range 
of angles considered, several diffraction minima are explored while the size 
of the distribution falls down by many orders of magnitude. The dotted 
line represents the result with no final-state interactions (PWIA), which is 
of course identical in both methods. The solid and the dashed lines are 
the results of method I and II, respectively. The two angular distributions 
are rather well correlated in all the kinematics here explored, except in the 
diffraction minima for p = q = 1.4 GeV/c. We have checked that these 
discrepancies are smoothed with increasing p, q until they almost disappear 
at p = q = A GeV/c, as it is clear in fig. 4. In any case, through all the 
kinematics considered the oscillatory patterns are very close to each other 
across a remarkably large range of variation in size. 

It must be noticed that the rich diffractive pattern of the angular distri- 
butions is partially due to the nontrivial structure of the PWIA contribution, 
which itself contains many local minima. We have already shown in a pre- 



vious work |Tj] that if the Woods-Saxon bound state \E'o is substituted by a 
pure harmonic oscillator, so to produce an exponentially decreasing angular 
distribution in PWIA, the results of both methods I and II still show an 
oscillatory pattern at large angles due to FSI. Thus, the natural interpre- 
tation is that the diffractive minima, which are reminiscent of the angular 
distribution for elastic proton-nucleus scattering , derive from the fact that 
the ejected proton is testing coherently the residual nucleus. This is pecu- 
liar of a completely exclusive reaction, where the residual nucleus does not 
fragment. Energy-integrated distributions (i.e. for a semi-inclusive (e,e'p) re- 



action pO|, ^) are by definition less sensitive to the structure of the recoiling 



{A — 1) system, thus leading to very different angular shapes. 

We have shown that for the ^^C(e,e'p)^^Bsi/2 reaction and for outgoing 
proton momenta in the range 1 < p < 4 GeV/c (relevant to the planned 
experiments at CEBAF) the eikonal approximation to the scattering wave of 
the ejectile produces FSI effects very similar to the ones obtained when the 
complete second-order differential equation is solved exactly up to 120 partial 
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waves. The angular distributions are in good agreement up to very large an- 
gles, where the absolute size can fall down by many orders of magnitude. The 
observed oscillatory pattern can be interpreted as a coherent diffractive scat- 
tering between the ejectile and the residual nucleus. Therefore, completely 
exclusive (e,e'p) reactions are best suited to verify this prediction. 

We would like to thank O. Benhar, S. Boffi, S. Jeschonnek, N.N. Nikolaev 
and S. Simula for many stimulating discussions. 
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Captions 



Fig. 1 - The distorted spectral density >S'|^/2 ^ ^ function of the depth U of 
the real part of the optical potential for the ^^C(e,e'p)^^Bsi/2 reaction 
at the outgoing proton momentum p — lA GeV/c and at the missing 
momentum pm = 0. The dashed line corresponds to the depth W = Q 
of the imaginary potential, the solid line to = 50 MeV, the dotted 
hne toW = 100 MeV. 



Fig. 2 - The distorted spectral density <S'|^/2 ^ ^ function oi p = q for the 
^^C(e,e'p)^^Bsi/2 reaction at pm = and with the optical potential 
depths [/ = 0, = p 50/1400 MeV. 



Fig. 3 - The distorted spectral density 5'|^/2 ^ ^ function of the angle 7 between 
p and q for the ^^C(e,c'p)^^Bsi/2 reaction with p = q = 1.4 GeV/c. The 
dotted line shows the PWIA result. The solid and dashed lines are the 
results of methods I and II, respectively (see text). 



Fig. 4 - The same as in fig. 3, but for p = g = 4 GeV/c. 
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